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I.  INTRODUCTION 


The  105mm  XM797  is  a  proposed  training  round  for  the  105mm  M735  high  veloc¬ 
ity  projectile.  The  requirements  for  this  training  round  pertinent  to  this 
study  are:  (1)  closely  simulate  the  trajectory  of  the  M735  for  ranges  up  to  3 
km  and  (2)  not  exceed  a  maximum  range  of  8  km.  In  order  to  achieve  this  per¬ 
formance,  the  XM797  was  designed  with  a  nose  cap  made  of  zinc  alloy  which,  due 
to  structural  weakening  caused  by  aerodynamic  heating  and  centrifugal  loads 
due  to  spin,  is  supposed  to  fail  and  result  in  the  projectile  breaking  apart. 
The  projectile  main  body  is  segmented  into  three  parts  which  are  held  together 
by  the  zinc  nose  cap.  A  simplified  drawing  illustrating  the  projectile 

concept  and  functioning  is  shown  in  Figure  1. 

Firing  tests  were  conducted  at  Yuma  Proving  Ground,  Arizona,  and  Fort 
Greely,  Alaska,  in  which  expected  performance  goals  were  not  met.  In  order  to 
obtain  a  better  understanding  of  the  behavior  of  the  zinc  alloy  nose  cap  for 
the  test  conditions  encountered,  a  computational  study  was  initiated  to  inves¬ 
tigate  the  in-depth  thermal  response  to  aerodynamic  heating  for  proposed  nose 
cap  configurations.  The  computational  technique  used  is  the  ABRES  Shape 
Change  Code  (ASCC).1’2’3  This  report  summarizes  the  results  of  a  computation¬ 
al  study  of  the  effects  of  aerodynamic  heating  for  several  zinc  alloy  and 
steel  nose  cap  configurations  pertinent  to  the  flight  conditions  for  the 
development  testing  of  the  XM797  shell. 


A.  Overview 


II.  COMPUTATIONAL  TECHNIQUE 


The  ABRES  Shape  Change  Code  (ASCC)  has  been  developed  over  a  number  of 
years  at  a  cost  of  tens  of  thousands  of  man-hours  and  several  millions  of 
dollars.  The  primary  purpose  of  the  ASCC  is  to  numerically  model  the  shape 


1.  Dahn,  T.J. ,  Cooper,  L.,  Rafinejad ,  D.,  Youngblood,  S.B.,  and  Kelly,  J.T., 
"Passive  Nosetip  Technology  (PANT  II)  Program.  Volume  I.  Inviscid  Flow 
and  Heat  Transfer  Modeling  for  Reentry  Vehicle  Nosetips,"  SAMSO-TR-77-11 , 
Space  and  Missile  Systems  Organization,  Air  Force  Systems  Command,  Los 
Angeles,  California  90009,  October  1976. 

2.  Rafinejad,  D.,  Dahn,  T.J.,  Brink,  D.F.,  Abbett,  J.J.,  and  Wolf,  C.J., 

" Passive  Nosetip  Technology  (PANT  II)  Program.  Volume  II.  Computer  Users 
Manual:  ABRES  Shape  Change  Code  (ASCC),"  SAMSO-TR-77-11 ,  Space  and 
Missile  System  Organization,  Air  Force  Systems  Command,  Los  Angeles, 
California  90009,  October  1976. 

ii.  Sandhu,  S.S.,  and  Murray,  A.L.,  "Reentry  Vehicle  Technology  (REV-TECH ) 
Program.  Volume  LIT.  Improved  Capabilities  of  the  ABRES  Shape  Change 
Code  (ASCC  79),"  Acurex  Report  TR-79-1 0/AS,  Acurex  Corpoixition/Acy'otherm , 
'186  Clyde  Avenue,  Mountain  View,  California  94042,  prepared  for  Sjxiee  and 
Missile  Systems  Organization,  Air  Force  Systems  Command,  Los  Angeles, 
California  90009,  duly  1979. 
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history  am.)  in-depth  unsteady  temperature  response  of  an  ablating  reentry 
vehicle.  The  code  includes  inviscid  and  boundary  layer  models  to  calculate 
heat  transfer  distributions  and  recession  rates  for  a  variety  of  flow  condi¬ 
tions. 

A  simplified  flow  diagram  for  the  ASCC  is  shown  in  Figure  2.  BRL  has 
sponsored  two  studies*4’5  which  were  performed  by  Acurex/Aerotherm  Corporation 
to  apply  the  ASCC  to  predict  the  aerodynamic  heating  response  of  high  velocity 
Army  shell  and  to  provide  a  basis  for  further  studies  to  be  carried  out  in- 
house  at  BRL  and  LCWSL. 

R.  Heat  Conduction 


The  heat  conduction  solution  is  carried  out  in  two  grid  systems:  (1)  a 
moving  body  oriented  grid  (s,r)  which  is  used  for  the  ablating  outer  layer  of 
the  body  and  (2)  a  fixed  cylindrical  coordinate  system  for  the  remainder  of 
the  body.  A  simplified  computational  grid  for  the  unsteady  heat  conduction 
computations  is  shown  in  Figure  3.  The  computational  technique  allows  each 
node  to  be  flagged  for  ident i f icat ion  of  material  properties.  This  allows 
multiple  materials  to  be  accounteu  for  in  a  simple  manner.  The  outer,  moving 
layer  is  solved  for  using  an  implicit  finite-difference  scheme.  The  inner 
region  is  solved  for  using  an  explicit  finite-difference  scheme.  The  equation 
solved  in  the  outer,  implicit  layer  is  Equation  (1). 


r  3T 

pCp  It 


1 _  |  3  r  /  b  \  dT  -l 

r„  (1  t('Ur/rc>  JSJ 


(1) 


d  r  /»  r  \  3 T -i  i  r  9  d i 

3r  rb  r  K  3r^  ^  P  p  n  3r 


3T 


where 


Cp  =  specific  heat 


body  circumferential  radius  of  curvature 


r0  +  r  cosO 


4.  Suehsland,  K.E.,  "Aerothernal  Aesessnent  of  Projectiles  Using  the  APSES 
Shape  Change  Code  ( ASCC) Acurex  Report  TH-80-31-AS ,  dune  1980. 

b .  Abbott,  M.J.j  Duiven ,  R.P. ,  Taub,  3.,  and  Beak,  R.A.S.,  "Thenml  and 

Structural  Analysis  of  draining  Round  Rose  Caps } "  Cont.mot.  Report  ARBRL- 
CR-004SL,  U.S.  Arrrry  Ballistic  Research  Laboratory /ARRADCOM,  Aberdeen 
Proving  Ground,  Maryland  2100b,  May  1981  (AD  A10071S). 


rc  -  local  streamwise  radius  of  curvature 
*  -  thermal  conductivity 

p  =  density 

n  -  surface  normal  recession  rate,  n  =  -r 
T  -  temperature 

t  =  time 

6  -  angle  between  normal  to  local  surface  and  axis  of 

Symmetry 

s  -  streamwise  distance  along  body 

r  =  distance  normal  to  body  surface  at  s 

The  equation  solved  in  the  inner,  explicit  region  is  Equation  (2). 


pC 


aT_ 

p  at 


3x 


i  3T .  1  3  .  oT , 

(K  +  y  37  (yk  37) 


(2) 


In  this  formulation,  the  lateral  heat  flux  is  assumed  to  be  small  compared  to 
the  heat  flux  in  the  normal  direction.  The  temperatures  of  computational 
nodes  on  the  same  ray  are  linked  implicitly.  The  temperatures  on  adjacent 
rays  are  linked  explicitly.  A  very  complex  iterative  procedure  is  used  to 
solve  for  the  surface  temperature  and  the  recession  rate  normal  to  the  outer 

surface,  n.  This  procedure  accomplishes  the  coupling  of  the  in-depth  solution 
with  surface  reactions,  convection,  and  radiation. 

C .  Inviscid  Flow 

The  inviscid  flow  field  provides  input  for  computation  of  the  outer  bow 

shock  shape  and  the  boundary  layer  development.  The  inviscid  flow  field  is 
modeled  in  ASCC  in  three  parts:  (1)  surface  pressure  distribution,  (2)  shock 

shape  and  shock  standoff  distance,  and  (3)  flow  field  details  in  the  shock 

layer.  Surface  pressure  distributions  establish  the  boundary  layer  edge  veloc¬ 
ities  while  pressure  levels  are  used  in  the  aetermination  of  thermodynamic 
and  transport  properties.  The  shock  shape  affects  the  heat  and  mass  transfer 
rates  through  curved  shock  effects  on  the  boundary  layer  edge  conditions. 

The  ASCC  surface  pressure  distributions  are  calculated  using  semi- 
empirical  correlations  established  for  supersonic  and  hypersonic  Might  of  re¬ 
entry  vehicle  nosetip  configurations.  The  correlations  are  evaluated  in  pach 
of  three  regions  (A.K.C)  of  the  flow  field  as  indicated  in  Figure  4.  These 

separate  results  are  then  combined  to  yield  the  complete  pressure  distrifcu 
tion. 


9 


The  location  of  the  sonic  point  Is  calculated  as  the  first  step.  After 
the  sonic  point  Is  located,  the  correlations  are  used  to  calculate  the  surface 
pressure  distribution  in  each  of  the  three  regions.  This  involves  application 
of  modified  Newtonian  relations  and  curve  fits  of  finite-difference  solutions 
for  inviscid  flow.  The  ASCC  code  also  allows  the  user  to  input  tables  of 
previously  determined  surface  pressure  distribution  if  desired. 

The  shock  shape  is  predicted  based  on  thin  shock  layer  assumptions  using 
an  integral  method  adapted  from  the  method  of  Crowell?  This  method  uses 
decoupled  integral  forms  of  the  continuity  and  axial  momentum  equations  along 
with  the  specified  surface  pressure  distribution.  Using  oblique  shock  rela¬ 
tions  and  the  equation  of  state,  conditions  behind  the  shock  are  evaluated  as 
functions  of  known  free- stream  conditions.  The  local  shock  angle  and  the 
shock  standoff  distance  are  the  two  unknowns.  The  resulting  nonlinear  equa¬ 
tions  are  solved  using  an  interative  technique.  The  relatively  small  nose 
radii  and  low  Mach  (Mach  <  5)  numbers  of  interest  to  the  Army  are  outside  the 
range  of  applicability  of  thin  shock  layer  theory  and  the  ASCC  method  can 
yield  shock  angles  and  shock  coordinate  values  which  are  inconsistent.  These 
inconsistencies  were  great  enough  to  cabse  the  code  to  fail  to  reach  a  solu¬ 
tion  for  some  geometries. 

A  modification  was  made  which  allows  the  shock  shape  to  be  predicted  by 
quadrature  of  the  shock  angle  alone.  This  procedure  gives  consistent  shock 
angles  and  shock  coordinates  for  the  shapes  considered  in  this  study. 

Another  difficulty  encountered  with  the  ASCC  shock  layer  modeling  was  the 
technique  for  determining  where  the  outward  normal  from  the  body  surface 
intersected  the  bow  shock.  For  slender  shapes  such  as  the  XM797,  the  original 
coding  could  result  in  "kinked"  shock  waves  that  were  multivalued  both  axially 
and  radially.  An  improved  calculational  procedure  was  developed  which  over¬ 
comes  this  difficulty  and  has  been  incorporated  by  Acurex/Aerotherm  in  updated 
versions  of  ASCC. 

D.  Viscous  Boundary  Layer  Modeling 

The  viscous  boundary  layer  determines  the  wall  heat  and  mass  transfer  and 
the  wall  shear.  The  predictive  method  employed  is  sufficiently  sophisticated 
to  include  all  the  important  physical  events  present  in  the  flow  field  (as 
long  as  there  are  no  embedded  shock  waves  in  the  flow  field)  and  is  computa¬ 
tionally  fast.  The  basic  concept  was  formulated  by  Dahn  et  al 1  and  developed 
Into  the  momentum-energy  integral  technique  (MEIT).  The  method  retains  the 
essential  features  of  the  boundary  layer  equations  and  accounts  for  physical 
events  such  as  blowing,  acceleration,  roughness,  property  variation,  and 
inviscid  vortical  layer  effects.  The  boundary  layer  methodology  is  the  most 
Important  individual  feature  of  ASCC.  The  basic  governing  equations  and  the 
methodology  are  outlined  below. 


0.  Cmmll,  P.G.,  "Shock  Shape  Calculation  in  Hody  Oriented  Coordinatin'  for 
Specified  Surface  Prcnourc, "  Acronpaco  Corporation  Memo  76-h  b-Olfl,  Mau 
Vr>n. 
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The  boundary 


layer  momentum  and  energy  integral  equations  are 


dp 


(I/O  h  <rpeue2Q)  -  Tw  +  (pv)wUi,w  +  "9?i 


(3) 


and 


U/r)  ~  (rVJe(ht,e  *  \U)  =  Qw  +  (pv)w(V  ,w  ‘  hw}  * 


(M 


The  various  quantities  that  appear  in  the  above  integral  equations  are 

defined  as: 


momentum  thickness  0: 


u  - 


u  -  u 


dy 


(3) 


energy  thickness  $: 


-  lv 


dy 


(6) 


boundary  layer  shape  factor: 


H  =  6*/e 


(7) 


displacement  thickness  6*: 


6*  - 


/’ 


(1  - 


pu 

P  u  • 
y)  1 


)  dy 


/  r>  \ 

(3) 


The  subscript  e  refers  to  properties  evaluated  at  the  boundary  layer  edge;  the 
ri  f  ‘  refers  !  in.iscid  properties  evaluate.,  along  the  stagnation  pe  nt 
T^nirops.  I  he  guantit,  u,,.  IS  the  inviscid  veloclt,  at  the  -all.  Phe 

terms  h.  and  ht  e  refer  to  local  and  edge  Inetscid  stagnation  entha.pies, 

respective.,,  »hich  ace  the  same  for  atmospheric  flight  )  . 


tT  JPV  w  U  ■  »  -  J  J  . .  -  - 

inviscid  outer  flow,  u,  =  ue,  f't,i  =  ht>e  ana  Pi  -  Pe 
equations  reduce  to  standard  integral  boundary  layer  equations. 
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To  solve  the  equations,  the  properties  at  the  edge  of  the  boundary  layer 
are  needed.  This  input  is  obtained  by  using  the  invlscid  flow  solutions  and 
an  entrainment  relation  which  accounts  for  the  entropy  change  due  to  shock 
curvature.  Dahn  et  al,1  based  on  theory  and  experimental  data,  developed  cor¬ 
relations  of  skin  friction,  Stanton  number,  and  auxiliary  relations  for  solv¬ 
ing  the  governing  equations.  The  auxiliary  relations  consist  of  a  series  of 
influence  coefficients  which  account  for  complex  Interaction  effects  between 
physical  processes.  The  Influences  modeled  include  Reynolds  number,  Mach 
number,  boundary  layer  conditions  (laminar,  transitional,  or  turbulent),  sur¬ 
face  roughness,  mass  transfer,  and  flow  field  acceleration.  The  solution  pro¬ 
cedure  consists  of  two  steps:  (1)  Series  solutions  at  and  in  the  vicinity  of 
the  stagnation  point  are  obtained.  (2)  A  finite-difference  scheme  is  used 
to  obtain  solutions  away  from  the  stagnation  region.  Initial  information 
required  includes:  (1)  surface  shape,  (2)  surface  pressure,  (3)  blowing  rate, 
and  (4)  temperature  distribution. 


A.  Overview 


III.  COMPUTATIONAL  RESULTS 


The  results  presented  here  represent  a  small  portion  of  the  total  computa¬ 
tional  effort  expended  in  this  study  of  the  aerodynamic  heating  for  the  XM797 
projectile.  The  intent  here  is  to  summarize  the  results  in  a  manner  that  will 
preserve  the  overall  scope  of  the  results,  indicate  the  modeling  capability 
inherent  in  the  ASCC  code,  and  discuss  uncertainties  discovered  in  comparison 
of  the  computational  and  experimental  results. 

It  may  be  surmised  that  the  input  data  to  ASCC  is  quite  complex  and 
requires  substantial  background  and  experience  to  select  correct  values  for 
the  many  input  parameters.  In  order  to  short  circuit  the  time  required  to 
gain  this  experience,  a  STAS  (Short  Term  Analysis  Services)  contract  was  car¬ 
ried  out  with  Acurex5  to  exercise  the  code  for  several  Army  shell  applica¬ 
tions.  This  served  as  the  starting  point  for  the  computational  studies  car¬ 
ried  out  at  BRL  and  LCWSL  for  the  XM797  projectile. 

The  results  will  cover  basically  three  configurations  of  the  XM797  nose 
cap:  (1)  all  zinc,  (2)  plug  nose  (steel  tip  with  zinc  afterbody),  and  (3)  all 
steel. 

B.  All-Zinc  Nose  —  Thermal  Simulation  Tests 


The  all -zinc  nose  represents  an  interesting  challenge  for  the  ASCC  since 
the  ability  to  model  ablation  with  the  accompanying  shape  change  is  required. 
An  example  of  the  grid  system  used  for  the  XM797  is  shown  in  Figure  5.  The 
locking  pin  protrusion  into  the  nose  cap  is  modeled  by  flagging  appropriate 
grid  nodes  for  the  material  properties  of  mild  steel.  The  location  of 
” critical  node  B"  is  indicated.  This  node,  was  selected  because  it  was  judged 
to  be  in  the  area  of  expected  failure  for  the  nose  cap.  (Critical  nodes  "A" 
and  "C"  are  discussed  later.)  Although  the  actual  failure  position  is  not 
known,  this  node  served  as  a  basis  for  comparison  between  computational 
results  and  experimental  test  data.  A  series  of  computations  was  carried  out 
to  examine  the  influence  of  the  computational  grid  resolution  on  the  predicted 
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in-depth  temperature  response  and  ablation  of  the  nose  cap.  The  results 
indicated  a  decided  sensitivity  for  surface  temperature;  however,  the  in-depth 
temperature  at  the  critical  node  was  not  strongly  affected.  The  results  to  be 
shown  here  used  the  maximum  array  capacity  built  into  the  code.  The  implicit 
grid  was  10  x  29  and  the  explicit  grid  was  60  x  25. 

Melting  produces  a  layer  of  liquid  at  the  surface  whose  flow  along  the 
body  is  not  considered  in  the  ablation  modeling  of  ASCC.  Although  modeling  of 
the  actual  mass  distribution  of  the  melt  was  not  attempted,  it  was  possible  to 
extend  empirically  the  existing  representation  of  transition,  turbulent  sur¬ 
face  roughness,  and  laminar  heating  augmentation.  This  was  done  by  specify¬ 
ing:  (1)  a  melting  temperature  for  each  material,  (2)  temperatures  at  which 
melting  induces  transition,  (3)  the  temperature  at  which  laminar  heating  is 
augmented  by  a  specified  amount  due  to  melting,  and  (4)  the  temperature  at 
which  turbulent  roughness  heating  *■:  augmented.  All  of  these  phenomena  are 
assumed  to  be  caused  by  melting  but  their  corresponding  temperatures  need  not 
be  the  melt  temperature.  These  degrees  of  freedom  allow  additional  adjustment 
of  the  model  to  improve  agreement  with  experimental  data  on  melting  surfaces. 

t 

After  analyzing  two  data  sources  on  zinc  noses  of  different  designs  (NASA 
Ames  thermal  wind  tunnel  tests7  and  actual  firing  tests8),  it  was  determined 
that  the  recession  rate  of  an  all-zinc  nose  was  closely  matched  by  an  effec¬ 
tive  turbulent  roughness  of  molten  zinc  of  about  20  mils,  a  laminar  heating 
augmentation  factor  of  10,  and  fixing  transition  at  the  most  forward  point  on 
the  zinc  surface  that  had  melted.  In  practice,  it  was  found  necessary  for 
numerical  stability  to  apply  a  ramp  in  roughness  and  heating  augmentation 
rather  than  a  step  at  the  specified  temperature  for  the  physical  phenomenon. 
This  ramp  extends  over  roughly  5°  F  on  both  sides  of  the  specified  tempera¬ 
ture. 

Some  experiments  have  been  conducted  in  the  AVCO  arc  jet  on  all-zinc 
noses.9  Figure  6  shows  the  stagnation  point  recession  as  a  function  of  time 
for  the  runs  available  (12357  and  12370).  Repeatabi  1  i ty  of  the  test  results 
was  not  good.  The  variation  may  be  due  to  the  random  size  and  location  of 
voids  in  the  cast  zinc  noses.  There  are  X-ray  data  and  sectioned  noses  that 
confirm  the  existence  and  prevalence  of  such  voids. 

Computational  predictions  using  the  previously  determined  20  mils  turbu¬ 
lent  roughness  follows  the  early  trend  of  Run  12357  well,  particularly  the 
time  of  onset  of  nose  recession,  but  falls  below  the  experimental  data  after 
one  second.  An  equivalent  roughness  at  melt  of  50  mils  was  also  tried  (see 


7.  Schwind,  Richard  G.,  "Hypersonic  Wind  Tunnel  Tests  of  Hose  Cap  Models 
Utilizing  Shape  Change  for  Range  Control,  "  Nielsen  Engineering  and 
Research ,  Inc.,  NEAR  TR  184,  March  1979. 

A' 

8.  Fleming,.  G.C.,  "Maximum  Range  Test  Results  of  the  105mm,  XM797  Training 
Projectile  Conducted  Between  January  and  October  1980  at  Yuma  Proving 
Ground,  Arizona ASD  IR-2-82,  December  1981. 

9.  Loeb,  A.,  LCWSL  unpublished  experimental  data.  Private  communication. 
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figure  6)  and  comes  closer  but  Is  still  below  the  experimental  results.  Run 
12370  has  the  same  time  of  onset  of  recession  as  Run  12357  hut  has  a  consis¬ 
tently  higher  recession  rate  and  is,  of  course,  even  further  above  the  predic¬ 
ted  recession. 

Since  more  detail  can  be  computed  in  the  nose  region  when  the  origin  of 
the  rays  of  the  polar  grid  system  is  closer  tu  the  nose,  the  origin  was  moved 
to  OX  -  0.6;  see  Figure  5  (using  50  mils  roughness).  The  computed  result, 
shown  in  Figure  6,  shows  recession  beginning  too  early  but  matching  the 
results  at  later  times  for  Run  12357  very  well.  Numerical  difficulties  are 
known  to  exist  with  the  current  grid  system  for  the  nose  reqion  of  a  slender 
body,  and  an  improved  system  is  currently  being  developed  on  contract.  The 
porosity  of  the  zinc  nose  is  another  unknown  and  unmodeled  effect  that  contrib¬ 
utes  to  the  discrepancies  described  above. 

Melting  and  ablation  are  similar  in  that  the  energy  balance  involves 
latent  heats.  However,  there  is  a  large  difference  in  the  magnitude  of  the 
latent  heats  for  melting  and  ablation.  Roth  ablation  and  melting  are  in-depth 
phenomena,  but  they  are  modeled  in  the  current  version  of  the  AS DC  code  as 
surface  phenomena.  This  difficulty  is  also  being  addressed  in  the  current 
contractual  effort. 

C .  A1 1  -7.i  nc  Nose  —  Atmospheric  Flight  Conditions 

Different  thermal /structural  analyses  had  resulted  in  different  conjec¬ 
tures  for  the  failure  points  and  modes  in  the  nose  (see  Figure  7),  the  lead¬ 
ing  candidates  being  shearing  or  opening  of  the  threads  in  the  zinc  at  point 
A,  bending  failure  at.  point  B,  or  hoop  stress  failure  at  point  C.  Therefore, 
an  attempt  was  made  to  see  if  a  unique  temperature  at  one  "critical  point" 
could  correlate  the  flight  tests  results  for  a  range  of  wall  thicknesses, 
launch  velocities,  initial  temperatures,  and  ambient  temperatures. 

The  velocity-time  trajectory  for  the  XM797  launched  at  standard  atmospher¬ 
ic  conditions  is  plotted  in  Figure  8.  An  example  of  the  predicted  shape  of 
the  nose  cap  at  several  times  in  the  flight  is  shown  in  Figure  9.  The  tem¬ 
perature  of  the  surface  at  one  second  of  flight  is  given  in  Figure  10.  The 
extent  of  the  melt  layer  is  clearly  indicated  by  the  region  of  constant 
temperature  of  1179  R. 

Examination  of  typical  computed  temperature  histories  the  points  of 
interest  (see  Figure  11)  shows  that  the  histories  are  quite  nonlinear  as  a 
function  of  tune.  An  exponential  rise  is  a  simple  form  that  approximates  the 
trend  of  the  actual  behavior. 


-  T  bt 
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The  form  ot  {.quit ion  (9)  was  used  to  correlate  predicted  temperatures  with 
observed  breakup  ^ failure)  times  for  points  A,  B,  and  C  for  selected  test 
data.  I  h  group  resulted  in  a  different  temperature  ot  breakup  at  the  same 
location,  with  maximum  spreads  of  20°  to  80°  F.  The  correlation  coef ( ii.ients 


14 


produced  by  Equation  (9)  for  the  3  locations  examined  were:  p, (A)  =  -0.970, 

Pj (B)  =  -0.467,  pj (C)  =  -0.963. 

The  longest,  experimental  tine  to  breakup  correlated  poorly  at  points  A, 
B,  and  C.  One  hypothesis  to  explain  this  result  is  that  molten  zinc  flows 
onto  the  cooler  steel  body  and  re-sol  idi  fies.  This  possibility  is  not  consid¬ 
ered  in  the  computational  modeling.  The  results  from  the  plasma  arc  jet 
tests  show  that  this  phenomenon  does  occur--model  disassembly  after  testing 
revealed  zinc  solidified  inside  the  split  aft  body  segments. 

The  important  correlation  is  predicted  breakup  time  versus  experimental 
breakup  time.  The  mean  temperature  at  breakup  was  used  to  obtain  the  time  of 
predicted  breakup.  These  results  are  plotted  in  Figures  12  to  14.  The  corre¬ 
lation  coefficients  for  these  cases  are:  Pt(A)  =  0.717,  Pt(B)  =  0.550,  pt(C) 

=  0.867.  Point  C  at  the  rear  interface  between  the  nose  and  the  body  clearly 
correlates  best.  The  mean  temperature  at  breakup  time  there  is  544K.  The  arc 
jet  tests  at  AVC0  clearly  showed  failure  at  point  C  in  the  hoop  stress  mode. 

0.  P 1  uy  Nose  Model  mg 

The  plug  nose  geometry  was  devised  in  an  effort  to  overcome  two  of  the 
difficulties  identified  in  testing  of  the  all-zinc  configuration:  (1)  higher 
drag  due  to  nosetip  ablation  and  (2)  inconsistent  performance.  A  simplified 
drawing  of  the  plug  nose  configuration  is  given  in  Figure  15.  The  shape  con¬ 
sists  of  a  steel  tip  and  an  afterbody  with  a  zinc  ring.  The  proposed  failure 
mode  still  relies  on  the  melting  or  thermal  weakening  of  the  zinc  due  to  aero¬ 
dynamic  heating. 

Computations  were  made  for  various  surface  rou, messes,  differing  boundary 
layer  transition  criteria,  and  trajectories  with  different  initial  projectile 
temperatures.  Examination  of  BRL  experimental  data  (aerodynamic  range  spark 
shadowgraphs)  indicated  that  the  ASCC  was  predicting  boundary  layer  transition 
too  early  cased  on  surface  length  and  free-stream  conditions.  The  experimen¬ 
tal  data  indicate  that  the  transition  Reynolds  number  (Re?)  should  be  about 

4.3  x  10b  whereas  the  code  predicts  Res  *  1.3  *  I0b.  The  code  also  allows  for 

prediction  of  transition  using  a  momentum  thickness  criterion.  Figure  16  shows 
an  example  of  the  sensitivity  of  the  in-depth  temperature  response  of  the  crit¬ 
ical  node  (B)  to  the  location  of  boundary  layer  transition.  Using  the  cri¬ 
teria  built  into  the  code,  the  critical  node  is  predicted  to  reach  melt  tem¬ 
perature  within  roughly  two  seconds  of  flight.  Forcing  transition  to  occur 
further  downstream  results  in  the  critical  node  never  reaching  the  melt  tem¬ 
perature.  These  computations  indicate  that  boundary  layer  transition  can  have 
a  significant  effect  on  the  in-depth  temperature  response  of  the  XM797  config¬ 
uration.  The  variation  in  predicted  heat  transfer  coefficient  at  time  one 
second  for  several  transition  criteria  is  shown  in  Figure  17.  The  location  of 
boundary  layer  transition  is  easily  identified  by  the  increase  in  heat  trans¬ 
fer  coefficient. 

A  series  of  experimental  firings  for  different  projectile  preconditioning 
temperatures  (cold:  244K,  standard:  294K,  and  hot:  322K)  was  performed  at 
Yuma  Proving  Ground,  Arizona.  Using  a  boundary  layer  transition  criterion  of 
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Res  =  1.75  ^  I0fc,  a  series  of  computational  results  was  obtained  for  flow 

conditions  matching  the  firing  tests.  A  comparison  between  the  computational 
results  for  the  critical  node  temperature  history  and  the  experimentally 
observed  projectile  breakup  times  Is  given  in  Figure  18.  The  notation  “B/U" 
indicates  time  of  projectile  breakup  as  determined  from  Doppler  radar  data. 
This  figure  shows  that  the  experimental  observations  correlate  qualitatively 
with  the  computational  predictions.  A  failure  criterion  is  implied  for  a 
critical  node  temperature  of  about  556K. 

E.  Solid  Steel  Nose  Modeling 

A  series  of  computations  was  run  for  nose  caps  made  of  solid  steel  in  an 
effort  to  identify  in-depth  regions  of  high  temperature  gradients.  Examples 
of  these  computational  results  are  shown  as  temperature  contours  in  Ficures 
19a  through  19f  for  specific  times  in  the  trajectory.  The  results  indicate 
that  very  high  temperature  gradients  occur  only  near  the  leading  edge  of  the 
projectile  and  that  the  temperature  gradients  become  noticeably  less  severe 
throughout  the  projectile  tip  after  one  second  of  flight. 

In  an  effort  to  evaluate  the  effect  of  the  location  of  boundary  layer 
transition  on  the  performance  of  the  XM797,  a  series  of  computations  was 
carried  out  for  different  locations  of  boundary  layer  transition  (X^).  The 

results  are  shown  in  Figure  20  where  the  time  for  the  critical  node  to  reach 
556k  is  plotted  versus  X-j-p  for  two  projectile  preconditioning  temperatures. 

These  results  indicate  an  inherent  sensitivity  of  the  design  concept  to  atmos¬ 
pheric  conditions  even  if  the  location  of  boundary  layer  transition  could  be 
fixed  using  a  boundary  layer  trip. 


IV.  SUMMARY 

This  paper  has  described  a  computational  study  in  which  the  ABRES  Shape 
Change  Code  (ASCC)  has  been  applied  to  predict  the  in-depth  temperature 
response  to  aerodynamic  heating  of  XM797  nose  cap  configurations.  Results 
have  been  shown  which  provide  examples  of  the  ability  of  the  code  to  predict 
effects  of  ablation,  location  of  boundary  layer  transition,  projectile  precon¬ 
ditioning  temperature,  arid  atmospheric  conditions.  Comparison  of  the  computa¬ 
tional  results  to  the  results  of  test  firings  provided  assistance  in  the 
analysis  of  the  projectile  performance,  including  identification  of  the  fail¬ 
ure  point  and  failure  mode  for  the  all-zinc  noses. 

The  experience  gained  in  the  course  of  this  study  has  resulted  in  the 
identification  of  several  desirable  modifications  to  the  code  which  would 


enhance  its  application  to  problems  of  interest  in  shell  design.  A  contract 
with  Acurex/Aerotherm  is  nearing  completion  which  will  result  in  the  code 
having:  (1)  an  improved  grid  configuration  for  long,  slender  shell;  (2)  an 

interactive  graphics  computational  grid  generation  routine;  (3)  improved 
modeling  of  in-depth  melt  and  contact  resistance  between  materials;  and  (4) 
capability  for  planar  2D  predictions  for  swept  fins. 


A  lesson  learned  in  this  study  is  that  aerodynamic  heating  cannot  reliably 
cause  a  functioning  late  in  the  flight  where  the  aerodynamic  heating  input  is 
reduced.  Any  desired  event  to  be  caused  by  aerodynamic  heating  must  occur 
early  in  the  flight  where  the  heat  transfer  rate  is  high. 


NOSE  CAP 


Figure  1.  Schematic  Illustration  of  the  105mm  XM797  Projectile  Concept 
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EXPLICIT  GRID 


Figure  3.  Simplified  Computat i ona  1  Grid  for  Heat  Conduction  Solution 


Figure  4.  Three  Regions  for  Inviscid  Flow  Modeling 
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Figure  7.  Simplified  Drawing  of  XM797  Nose  Cap  Showing  Critical  Points 
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17.  Distribution  of  Heat  Transfer  Coefficient  for  Several 
Criteria  for  Boundary  Layer  Transition 
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